function cons_sum_interfe = cal_sum_interfe(N,K,ww,hGA)
cons_sum_interfe=0;
for nn=1:N
    for kk=1:K
        cons_sum_interfe=cons_sum_interfe + (hGA{nn}'*(ww(:,nn,kk)*ww(:,nn,kk)')*hGA{nn});
    end
end
cons_sum_interfe = real(cons_sum_interfe);
end